Spiral charge frustration in molecular conductor (DI-DCNQI)2Ag 
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We theoretically study the effect of spiral-type charge frustration in a quasi-one-dimensional 
molecular conductor (DI-DCNQI)2Ag. We clarify how the spiral frustration in the interchain 
Coulomb repulsion is relieved and leads to a self-organization of complex charge-lattice ordered 
chains, in agreement with the recent synchrotron X-ray study [T. Kakiuchi et al, Phys. Rev. Lett. 
98, 066402 (2007)]. In addition, we find that a keen competition between charge and lattice degrees 
of freedom under the frustration gives rise to a characteristic temperature within the ordered phase, 
below which a drastic growth of molecular displacements occurs. Our results enlighten the relevance 
of the spiral frustration and provide a possible reconciliation among puzzhng experimental data. 

PACS numbers: 71.30.-|-h, 71.10.Fd, 71.20.Rv, 71.28.-|-d 
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Geometrical frustration in charge degree of free- 
dom, i.e., charge frustration, has been a fundamen- 
tal problem in condensed matter physics, since Ander- 
son argued its importance in the Verwey transition in 
magnetite [ij. Recently, it is attracting renewed in- 
terests partly because exotic behaviors are found in 
several molecular conductors, where frustrated lattice 
structures are often formed. A typical example is 
the bis(ethylenedithio)tetrathiafulvalene (BEDT-TTF) 
based compounds, in which strong frustration in the two- 
dimensional triangular lattice makes the charge ordering 
(CO) unstable [3| and intriguing transport phenomena 
are observed in the critical region d, Q . 

Quasi-one-dimensional (QID) molecular conductors, 
however, have been less examined from the viewpoint 
of frustration; it is usually supposed that the ID chains 
are well decoupled, offereing canonical systems to study 
ID physics such as Tomonaga-Luttinger liquids, charge 
and spin density- waves, and so on [5|. Nevertheless, 
the QID compounds also often retain a "hidden" three- 
dimensional (3D) frustration in packing the ID chains. 
Such interchain frustration effect has been studied within 
extended Hubbard-type models on frustrated zigzag lad- 
ders [i,0,i. 

We consider herea typical QID molecular con- 
ductor (DI-DCNQI)2Ag 0], which provides a dif- 
ferent avenue to interchain frustration. This com- 
pound consists of ID chains of DI-DCNQI (2,5-diiodo- 
iV, iV'-dicyanoquinonediimine) molecules whose elec- 
tronic structure near the Fermi level is described by a 
QID band at quarter filling [l^. A prominent feature 
of this system is the phase transition at Tc ~ 220 K. 
It was observed that the ^^C-NMR line splits below Tc, 
which was ascribed to a Wigner-crystal type CO with a 
twofold period along the chains Although the sim- 
ple CO scenario was supported by theoretical studies in 
ID models [l 2| . there have been many puzzling behaviors 
found in experiments. 

One is the lattice displacements suggested by optical 
spectra analyses [l^ [ij]: The IR and Raman spectra 
are complicated and not well fitted by such simple CO 



without assuming any lattice distortions. Another unac- 
countable behavior is seen in the NMR line shape 
There appears an extraordinarily wide broadening below 
^100 K in the outer line with a larger shift, which can 
hardly be ascribed to critical fluctuations related to the 
magnetic ordering at a considerably lower temperature 
(T) of 5 K. Furthermore, a small peak starts to grow at 
around zero shift in a similar T range. Such a character- 
istic T range within the ordered phase was also observed 
in transport properties. The 1/T derivative of the log 
of resistivity shows, in addition to a peak structure at 
Tc, a broad hump centered at around 100 K TB*]. In the 
same region, the dielectric permittivity exhibits a peak 
and large frequency dependence [l6|. All these behav- 
iors, implying intrinsic changes at a much lower T than 
Tc, are difficult to understand from the simple ID CO 
picture. 

A clue to a possible origin of these features was found 
by a recent experimental study based on synchrotron X- 
ray crystal structure analysis at 50 K 17]. It was claimed 
that the ordered structure is not a simple CO: Instead 
it is a periodic array of DI-DCNQI chains with mixed 
patterns of charge-lattice symmetry breaking as shown 
in Fig. [1] which we call the "mixed state." The impor- 
tance of frustration in the lattice structure was pointed 
out; if all the chains exhibit the simple twofold CO, they 
would suffer from a spiral-type charge frustration among 
neighboring chains, since the DI-DCNQI molecules are 
displaced with each other by 1 /4 of lattice spacing along 
the chain as shown in Fig.[T]Jb) (see also Fig. 4 of ref. ^17j). 
However, it is not clear how the mixed state is stabi- 
lized and, furthurmore, what causes the peculiar T de- 
pendences in the ordered phase mentioned above. 

Theoretically, interplay between charge and lattice de- 
grees of freedom in QID molecular conductors has been 
studied on the basis of ID or QID extended Hubbard 
models at quarter-filling [l^. In particular, recent stud- 
ies showed that CO and lattice dimerization (LD) com- 
pete with each other but they coexist in some cases [l^ . 
These three states, CO, LD, and their coexistence, are 
in fact realized but in different chains in the mixed state 
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FIG. 1: (color) The charge-lattice ordered state (the mixed 
state) in (DI-DCNQI)2Ag (a) DI-DCNQI and Ag sites 
viewed from the chain direction (c axis), where the arrange- 
ment of chains with different orders are indicated; CO and LD 
represent charge order and lattice dimerization, respectively, 
(b) Side view of the DI-DCNQI sites where the size of circles 
and arrows represent the charge density and the lattice dis- 
placement, respectively. The blue dotted lines in (b) are the 
spiral V' bonds in eq. ((T}. 

shown in Fig. [T] The spiral frustration is expected to 
be a key factor for the nontrivial mixture of orderings 
out of uniform and equivalent chains, however, theoret- 
ical studies taking account of it have not been done so 
far. 

In this study, we investigate the effect of the interchain 
spiral frustration in QID systems for a model of (DI- 
DCNQI)2Ag. We clarify how the mixed state in Fig. [T] 
is stabilized under the charge frustration. We also cal- 
culate finite-T properties of the charge-lattice coupled 
phenomena and propose a possible scenario to reconcile 
the puzzling experimental data in this compound. 

We consider a QID model which takes account of the 
3D interchain spiral structure in the actual compound. 
Each ID chain is represented by the extended Peierls- 
Hubbard model at quarter filling [l^. Here, in order 
to concentrate on the competition between charge and 
lattice degrees of freedom, we take the limit of strong 
onsite Coulomb interaction and leave out the spin degree 
of freedom [l^. Namely, we study a half-filled spinless 
fermion model coupled to the lattice, whose Hamiltonian 
reads 

where a\ (ui) is the creation (annihilation) operator 
for a spinless fermion at the ith DI-DCNQI site, and 
Ui = alui. Here, the first term is the ID part where 
the sum on (i,j)id is taken for nearest-neighbor sites 
along the chains. The next two terms are the inter- 
chain Coulomb interactions where are the nearest 
interchain site-pairs shown by the blue dotted lines in 




FIG. 2: (color) The CO-I-LD state. Notations are the same 
as in Fig.Ill whereas CO-fLD(l) and CO-(-LD(2) in (a) indi- 
cate the two kinds of chains, both with charge order as well 
as lattice dimerization but with different degrees. The thick 
V' bonds in (b) suffer from frustration in Coulomb repulsion 
between charge rich sites. 

Fig.[ljb), and (i,jy± are the next-nearest interchain pairs 
such as sites 1-8, 2-7, etc., while interactions of farther 
range are neglected for simplicity. Note that both terms 
undergo the spiral frustration. The transfer integrals and 
the Coulomb repulsions are coupled with the lattice dis- 
placement Ui of each molecule, as Uj = i(l -I- aSuij), 
= V{1 + f35u,j), Vlj = V'{1 + P'Su,,), and V^j = 
V"{1 + (3" Suij), where Suij — Ui — Uj [20]. The last term 
in eq. ([1]) is the elastic energy, and the displacements Ui 
are treated as classical variables within the adiabatic ap- 
proximation. In the following, we restrict ourselves to 
Ui being uniaxial, i.e., the molecules are allowed to move 
only in the chain direction, by considering the experi- 
mental results [13] ■ 

We treat the Coulomb interaction within Hartree-Fock 
approximation. The size of unit cell we consider is the ex- 
perimental one which contains 8 chains with 2 sites each, 
namely, 16 sites as shown in Fig. [TJ The order parame- 
ters, (ui), {a\aj)^ and Ui are self-consistently determined 
with the precision of 10~^. We set t = 1 and a — 1 (ab- 
sorbed in the notation of u^), and fix y = 1.5, V = 0.75, 
V" = 0.15, (3 = 0.5, /?' = 0.0325, and /?" = 0.0975 [IH, 
while varying K^^, which represents the softness of the 
lattice. We will comment on the choice for the intersite 
Coulomb repulsions later. 

First let us discuss the results for the ground state 
properties. We find that the charge-lattice ordering 
equivalent to the mixed state pattern in Fig. [1] is indeed 
stabilized over others in a certain range of K~^, among 
the several solutions self-consistently obtained. Other- 
wise, another solution gives the lowest energy: In this 
state, all chains show both CO and LD as shown in Fig. [21 
which we call here the CO-I-LD state. Both of these states 
are insulating. Figure [3] shows the ground state energies 
as well as the charge disproportionations 5i — {rii) — 1/2 
and the lattice displacements Ui for these two solutions 
as a function of K'^ . For K'^ < 0.10, the CO+LD 
state is the ground state: As shown in Fig.[3I^d), the lat- 
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FIG. 3: (color) The ground state properties of the mixed state 
and the CO+LD state. dependences of (a) the ground 

state energies, (b) energy difference between the two states, 
and (c), (d) charge disproportionations 5i — (rii) — 1/2 and 
lattice displacements Ui for representative sites in the two 
solutions. Si and Ui for other sites are equal or have opposite 
sign with equal absolute value to the plotted ones (see Figs. [T] 
and [2]). Note that (145,117) in (c) and (112,1*5) and (m7,wio) 
in (d) take almost the same but not identical values [see also 
the insets of Figs.|4jc) and (d)]. 

tice displacements Ui are almost negligibly small in this 
region; therefore, it is essentially regarded as a bundle of 
ID CO states stabilized by the intrachain Coulomb re- 
pulsion V. This state suffers from the charge frustration 
in the thick V' bonds in Fig. [HJb). On the other hand, 
for > 0.19, the CO+LD state is the ground state as 
well [l^l; when is large, the charge disproportiona- 
tions Si become small, and hence, it is essentially a LD 
state stabilized by the intrachain charge-lattice coupling. 
In between these two, the mixed state gives a lower en- 
ergy for 0.10 < K^^ < 0.19. This state emerges in the 
transient regime from the CO-dominant state for small 
to the LD-dominant state for large K~^, both of 
which undergo frustration. There, the system finds a 
way to relieve the frustration by self-organizing arrays of 
the CO, LD, and coexisting chains. 

Next we show how these coexistences and competitions 
evolve at finite T. In Fig.dl^a), T dependences of the free 
energies for the two solutions are shown at = 0.15, 
for which the ground state is the mixed state. As in- 
creasing T, the free energy difference between the two 
states becomes smaller, and at T^^^ = 0.66, a first-order 
transition takes place above which the CO-I-LD state 
gives lower free energy [see also Fig. Ul^b)]. Finally, for 
T > T^CO-i-LD ^ Q system turns into the uniform 

metallic state. The sequence of the phase transitions is 
summarized in Fig.OJa). 

An important point is that the free energy difference is 
extremely small for T^^*"* <T < T^CO-i-ld ^ j^^^h less than 
10~^ as shown in Fig. [UJb). Another key observation 
is that T dependences of Si and Ui are similar between 
these two solutions as plotted in Figs.^c) and (d). Their 
transition temperatures from the high-T uniform phase 
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FIG. 4: (color) Temperature dependences of (a) the free en- 
ergies, (b) free energy difference, and (c), (d) charge dispro- 
portionations Si — (n,;) — 1/2 and lattice displacements Ui for 
the two solutions for = 0.15. Insets show enlarged view 
for Ui near the transition temperatures. 

are almost the same (T 0.91), and by decreasing T, 
charge disproportionations Si develop first while lattice 
distortions Ui remain suppressed. When T is decreased 
down to T* ~ 0.7, then the Ui rapidly develop. Note that 
this crossover temperature T* almost coincides with T^^*. 
The first order phase transition as well as the finite-T fea- 
tures described above are always seen in our calculations 
when the mixed state is the ground state. 

From these results at finite T, we propose two differ- 
ent scenarios for what is taking place in the actual com- 
pound (DI-DCNQI)2Ag. One is that, as in our calcula- 
tion, the system exhibits two phase transitions. The first 
one is the transition at — 220 K, which corresponds 
to T = T^O+LD^ and the second occurs at T^^"* well be- 
low Tc as indicated in Fig. [S{a); yet this latter transi- 
tion has not been observed in experiments so far. The 
other scenario is that the CO-I-LD state in our calcula- 
tion is destabilized if fluctuations are taken into account 
beyond the mean-field level, and the mixed state gives a 
lower free energy in the entire T range within the ordered 
phase. Then, we predict that there is only one transition 
to the mixed state; Tc = 220 K corresponds to TJ"'™'^, 
which is followed by a crossover behavior at T*, instead 
of the successive transitions, as illustrated in Fig. \Slh). 
This latter scenario is based on the observation that CO 
suffering from frustration is often destabilized by fiuctu- 
ations [l^, , and the free energy difference is already 
very tiny in the present mean-field results. 

In both of these scenarios, an additional characteristic 

(a) _,ist CO+LD (b) jT-T * T- mixed 

.(c ^ T 

I — 1 1 ^ I — t1 1 >i 

mixed CO+LD mixed 

FIG. 5: (color) Schematic illustration of the proposed two 
scenarios for the phase transition(s) in (DI-DCNQI)2Ag. 
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temperature is predicted: T^^* in the former and T* in 
the latter. In either case, the system undergoes drastic 
changes there in the lattice sector, i.e., a rapid devel- 
opment of the lattice displacements within the ordered 
phase. This is suggestive in considering the puzzling ex- 
perimental results mentioned above. The rapid growth 
and the expected large fluctuations in the lattice displace- 
ments may explain a peak and frequency dependence of 
the dielectric permittivity [l6| as well as the anomalous 
broadening and an additional peak in the NMR spec- 
tra A hump in the 1/T derivative of the log of resis- 
tivity [l5| can be also related to the drastic changes at the 
characteristic temperature. To confirm our proposal, sys- 
tematic measurements of T dependence of X-ray diffrac- 
tion as well as NMR for single crystals are highly desired. 

Let us comment on the values of interchain Coulomb 
repulsions used in our calculations. By investigating sev- 
eral sets of parameters with varying K~^, we have found 
that the mixed state is robust particularly for V' ^ V", 
while the region becomes small when V' is compati- 
ble to V" . The results imply that, since the center- 
to-center distances between two DTDCNQI molecules 
are not much different for V and V" , other factors are 
also important in determining the (effective) interchain 
Coulomb repulsions, such as the cation environment and 
the anisotropy of molecular orbitals. In particular, the 
displacements of Ag+ cations to positions along the V' 
bonds can be relevant by enhancing the values of V 
(and reducing V") in our model, since the cation with 
positive charge enhances charge disproportionations be- 
tween the neighboring DI-DCNQI molecules. A similar 
effect was discussed in another QID molecular conductor 
(TMTTF)2X [24, 25]. 

Finally we make a remark on the spin degree of free- 
dom. Experimentally, an antiferromagnetic transition is 
observed by NMR [| and ESR at T ~ 5 K. In the 



mixed state, we expect localized spins to appear on the 
charge rich sites in the CO chains, on the dimers in the 
LD chains, and on the dimers in the coexisting chains but 
with more spin density on the charge rich sites. Both in- 
trachain and interchain exchange couplings depend on 
the degree of CO as well as LD, which show characteris- 
tic T-dependences as discussed above; hence, it is highly 
nontrivial how the spin degree of freedom acts under such 
condition. This should be a reason why simple Bonner- 
Fisher type analysis fails to understand the magnetic sus- 
ceptibility and specific heat data [13] ■ Calculations 
including spins are left for future study. 

In summary, we have studied a charge-frustrated spin- 
less fermion model coupled to the lattice degree of free- 
dom to understand the charge-lattice coupled phenom- 
ena in the quasi-one-dimensional molecular conductor 
(DI-DCNQI)2Ag. We have shown that a peculiar mixed 
state, where charge ordered, lattice dimerized, and their 
coexisting chains arrange periodically, is self-organized so 
as to relieve the spiral frustration between the chains. We 
have proposed two scenarios for the finite-temperature 
properties, and in both of them there is a characteris- 
tic temperature where lattice distortions develop rapidly. 
This drastic change provides a key to understand the puz- 
zles remaining in experiments. Similar interchain frustra- 
tion is seen in other quasi-one-dimensional compounds, 
and therefore, our results indicate the necessity to ex- 
amine its effect to obtain deeper understanding of these 
systems. 
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